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ABSTRACT 

We have performed the first measures of mass accretion rates in the core of the Orion Nebula Cluster. 
Four adjacent fields centered on the Trapezium stars have been imaged in the U- and B-bands using the 
' Wide Field Planetary Camera 2 on board the Hubble Space Telescope. We obtained photometry for 91 

O '■ stars in the U-band (F336W) and 71 stars in the B-band (F439W). The WFPC2 archive was also searched 

' to obtain complementary V-band (F547M) and I-band (F791W) photometry. In this paper we focus our 

(N ■ attention on a group of 40 stars with known spectral types and complete UBVI WFPC2 photometry. 

. We locate each star on the HR diagram considering both the standard ISM reddening law with Ry = 3.1 

D ' and the "anomalous" reddening law with Ry — 5.5 more appropriate for the Orion Nebula. Then we 

pL^ ', derive the stellar masses and ages by comparing with the evolutionary tracks and isochrones calculated 

by D Antona & Mazzitelli and Palla & Stabler. Approximately three quarters of the sources show excess 
luminosity in the U-band, that we attribute to mass accretion. The known correlation between the 
, U-band excess and the total accretion luminosity, recalibrated for our photometric system, allows us to 

^ ■ estimate the accretion rates, which are all found to be in the range 10~^ — 1O~^^M0 yr~^. For stars older 

lO i than 1 Myr there is some evidence of a relation between mass accretion rates and stellar age. Overall, 

0^ ■ mass accretion rates appear lower than those measured by other authors in the Orion flanking fields or 

' in Taurus- Auriga. Mass accretion rates remain low even in the vicinity of the 10~^ Mq yr~^ birth line 

of Palla & Stabler, suggesting that in the core of the Trapezium cluster disk accretion has been recently 
depressed by an external mechanism. We suggest that the UV radiation generated by the Trapezium 
OB stars, responsible for the disk evaporation, may also cause the drop of the mass accretion rate. In 
, this scenario, low-mass stars may terminate their pre-main sequence evolution with masses lower than 

' those they would have reached if disk accretion could have proceeded undisturbed until the final disk 

consumption. In OB associations the low-mass end of the Initial Mass Function may therefore be affected 
Q . by the rapid evolution of the most massive cluster's stars, causing a a surplus of "accretion aborted" very 

$H ' low-mass stars and brown dwarfs, and a deficit of intermediate mass stars. This trend is in agreement 

I with recent observations of the IMF in the Trapezium cluster. 

I Subject headings: open clusters and associations: individual (Orion Nebula Cluster) — accretion, 

^ . accretion disks — stars: pre-main-sequence — stars: formation 

•i-H . 

^ ' 1. INTRODUCTION 

One of the most relevant processes in the early stellar evolution is the interaction between young pre-main-sequence 
stars and their circumstellar disks. With the exception of the very initial phases of star formation, during which protostars 
may directly accrete low angular momentum material radially infalling from their parental cloud, most of the stellar mass 
build-up occurs through disk accretion (e.g. Tereby, Shu, & Cassen 1984). Disk accretion regulates the final stellar mass 
removing angular momentum via viscous dissipation and feeding conspicuous mass loss through collimated outflows (e.g. 
Hartmann 1998). There is evidence that the disk accretion process becomes less and less significant as the young star ages 
(Hartmann et al. 1998), and, for low mass stars {M < IMq), accretion must terminate within the disk lifetime (6-10 Myr, 
Haisch et al. 2001; Skrutskie et al. 1990), well before reaching the main sequence. However, a quantitative understanding 
of the evolution of the accretion process is still missing. In principle, mass accretion can be studied in a variety of ways. 
The infrared disk emission produced by internal viscous dissipation is potentially an important diagnostic tool, but in 
most disks direct irradiation of the surface from the central star is a stronger heat source (Kenyon & Hartmann 1987). 
A less ambiguous indicator is the recombination radiation produced when the disk material falls onto the stellar surface 
(Bertout 1989), through a boundary layer or, more probably, through accretion columns in the stellar magnetosphere 
(Konigl 1991; Camenzind 1990, Shu et al. 1994). The hot gas continuum adds up to the stellar continuum and alters the 
depths of photospheric lines, the so-called "veiling." From the amount by which stellar absorption features are veiled, one 
can make quantitative estimates of the relative emission of star and hot continuum, obtaining the integrated accretion 
luminosity and, therefore, the mass accretion rate (Edwards et al. 1987; Hartigan et al. 1991; Bathala et al. 1996; Hartigan, 
Edwards & Ghandour 1995; Gullbring et al. 1998). The hot continuum emission becomes easily detectable shortward of 
the Balmer discontinuity as a characteristic ultraviolet excess (Kuhi 1974). Gullbring et al. (1998) have shown that the 
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excess radiation in the Johnson U-band is related to the accretion luminosity derived from veiling measures. Calvet & 
GuUbring (1998) have reproduced the relation within the framework of the magnetically driven accretion column model, 
showing that the correlation does not depend on the main stellar parameters (mass, effective temperature) at least for 
spectral types in the range from M3 to K5. The UV excess has been used by Gullbring et al. (1998, 2000) to determine 
the mass accretion of stars in the Taurus- Auriga association, by Hartmann et al. (1998) in Chamaeleon-I, and by RebuU 
et al. (2000) on the largely unexplored outer regions of the Orion Nebula. 

This paper presents new HST photometry of the Trapezium cluster, the nearest region of ongoing high-mass star 
formation. Located in the core of the Orion Nebula, this cluster contains at least 3000 members between ~ 40 Mq 
and ~ 10 Mj (Lada & Lada, 2003). The availability of spectral type classification for « 1000 stars (Hillenbrand 1997) 
provides a unique opportunity to build the largest database of accretion luiniiiosities and mass accrcrtion rate. In this work 
we present the results obtained for a sample of 40 isolated point sources with known spectral types and accurate HST 
photometry. The observations are discussed in Section 2, whereas in Section 3 we illustrate the criteria leading to the 
sample selection. Our data analysis strategy is presented in Section 4. In Section 5 we present the results, in particular 
the location of our sources on the HR diagram and the derived mass accretion rates. Finally, in Section 6 we discuss our 
findings. Whereas our original goal was to build the first reliable map of the mass accretion rates vs. stellar mass and 
age, we find evidence of a link between the mass accretion rates and the extreme environmental conditions of the cluster 
core. This may have profound implications on the final stellar masses, affecting the Initial Mass Function (IMF) of the 
Cluster's core, especially at the low mass end. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. HST Data sets 

Our observations were made in March 2001 using the Wide Field Planetary Camera 2'* on board HST. Four fields 
centered in the immediate surroundings of the Trapezium stars were imaged with the F336W and F439W filters. The 
F336W filter has a long cutoff wavelength at A ~ 3600 A, shortward of the Balmer discontinuity, and is therefore similar 
to the Stromgren u filter. The F336W filter is affected by red-leak, and in Section 4.1 we discuss how this has been taken 
into account. The F439W filter closely matches the standard Johnson-B passband. For each field, we obtained four 400 
second exposures with the F336W filter and two 180 second exposures with the F439W filter. Short exposures of 1 and 
30 second were also taken to measure the bright stars. 

We complement these observations by using archival WFPC2 images of the Trapezium region taken with F547M and 
F791W broad-band filters, which are close to the standard Johnson V and Cousins I passbands. The largest number of 
counterparts to our sources is found in the data set of GO proposal 6666 (PI Stauffer). Ten fields were observed in November 
and December 1998, taking on each field two 300 second long exposures with the F791W filter, and two 500 second 
exposures in F547M, filter, complemented by a 60 second exposure per filter. Other data come from observations made 
in March 1995 for GO proposal 5469 (PI Bally), where four fields were imaged in the F547M filter, with nine 40 second 
exposures plus five 180 second exposure for one field, and three 30 second exposures on the remaining three. Finally, for 
GO proposal 6603 (PI Bally), four fields were observed in the F547M filter, with a total of three 60 second exposures 
and two 30 second exposures. The first HST observations of the Orion Nebula were made in 1991 (GO proposal 2595, PI 
Stauffer), using the Planetary Camera channel of the Wide Field Camera ( WFC). At that time the HST was still affected 
by spherical aberration. These data, published by Prosser et al. (1994), have been later included in main photometric 
catalog of the Orion cluster (Hillenbrand 1997), assembled assigning the highest priority to the observations made with 
the finest pixel sampling. We shall use them for those sources lacking more recent WFPC2 photometry. 

2.2. Data reduction 

The HST pipeline produces data with bias, dark current, and flat-field corrections applied. Using the stellar photometry 
program, HSTphot^, specifically designed for use with HST/WFPC2 images, we removed bad columns, cosmic-rays, and 
hot pixels. To avoid objects whose spectral type classification may be affected by non-stellar emission, we confine our 
study to single, well isolated stars, and exclude sources that appear close pairs, surrounded by photoionized envelopes, 
and in general resolved. It must be remarked, however, that even the most prominent dark silhouette disks in Orion can 
only be reveled in narrow band images, since in broad-band filters they are embedded under the wings of the Point Spread 
Function of the central star. There is abundant evidence from IR data that disks are present around at lest ~ 80% of 
the source in our fields (Hillenbrand & Carpenter 2000, Lada et al. 2000), and our sample is not biased against sources 
with circumstcUar disks. Astrometry and photometry have been performed on individual CCD images. We obtained 
coordinates for all stellar sources by using the IRAF/STSDAS task METRIC, which corrects for geometric distortion 
of WFPC2. According to the WFPC2 manual, the final relative positions are expected to be accurate to better than 
5 mas for targets contained within the same CCD chip, and O'.'l for targets on different chips. We used daophot for 
multi-aperture photometry and hstphot for psf photometry, obtaining consistent results. The hstphot data have been 
generally used for the analysis. The excellent quality of the residual images after PSF subtraction confirms that spurious 
contamination from the uniform nebular background has been virtually eliminated. The photometric zero points for each 

^ updated version of WFPC2 Instrument Handbook for Cycle 12 available at www.stsci.edu/instruments/wfpc2/Wfpc2_hand/wfpc2_handbook.html/ 
^ HSTphot manual available at www.noao.edu/staff/dolphin/hstphot/ 
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chip, filter, and gain setting were taken from Dolphin,^ relative to an aperture size of 0'.'5 (i.e. 5 pixels for the Wide Field 
frames and 11 pixels for the Planetary Camera frame) in radius. They are known with an accuracy of ~ 0.005™. 

3. THE SAMPLE 

On the basis of the criteria described in the previous section, wc identify 91 sources in the F336W images and 71 sources 
in the F439W images. For 49 sources we have both magnitudes, and for 40 of them we also have a spectral classification 
from Hillenbrand (1997). In this paper we concentrate our attention on these 40 objects. In Table 2 we present their 
photometry in the four WFPC2 filters. Since most of the stars with known spectral types are relatively bright, they 
appear saturated in the archival V- and I-band WFPC2 images. For this reason, only 26 stars have WFPC2 photometry 
in all four UBVI-bands. For the remaining objects with incomplete or absent WFPC2 V- and I-band photometry we rely 
on the catalog of Hillenbrand (1997). 

The relevance of stellar variability can be addressed using the multiple F547M archive images taken at different epochs. 
The last two columns of Table 1 report the number of F547M observations and the scatter of their values. In Figure 1 we 
plot the individual observations, showing for each star the scatter of all available measures vs. their mean value. The cross 
symbols represent the Hillenbrand (1997) V-band magnitudes, that for the brightest sources are the only available data. 
For 36 stars we have two, or more, measures, and the standard deviations range from a few hundredth of magnitudes to 
more than half a magnitude. Weighting each standard deviation by the corresponding number of measures, we obtain 
an average ct(F547M)= 0.211™. In 2/3 of the cases (21 stars) the peak-to-pcak variability is larger than 0.2™. Hcrbst 
et al. (2002) have recently presented the results of their study of stellar variability in the Orion Nebula cluster, done in 
a filter close to our F791W passband. Their study indicates that essentially all stars brighter than Ic 16™, i.e. with 
photometric accuracy better than 0.01™, are variable. In particular, about half of the stars show peak-to-peak variations 
larger than ~ 0.2™, in excellent agreement with our ratio. Variability directly affects the derived stellar luminosity and 
V — I color excess, sincx; the niciasures in the F791W and F547M filters arc in general not simultaneous. The resulting 
uncertainty on the extinction affects both the total and accretion luminosities. 

The V-magnitudes presented in Table 1 are the averages of all available measures, including the data of Hillenbrand 
(1997). The associated errors, on the other hand, are the average errors of the F547M measures made with WFPC2. The 
errors are therefore indicative of our measure accuracy and not of the uncertainty on the "true" stellar flux. Figure 2 and 3 
show the spectral energy distributions, together with the model results presented in the forthcoming sections. The near-IR 
photometry plotted for completeness has been taken from Hillenbrand et al. (1998) for the J-band, and Hillenbrand & 
Carpenter (2000) for the H- and K-bands. 

4. ANALYSIS 

To obtain the mass accretion rates, we adopt the following strategy. First, for each star we compare our V- and I- 
band observations with the reddened spectral energy distributions of template stars, deriving the extinction towards each 
individual star, and the bolometric luminosity. Then, we correct the observed flux in the U-band for the extinction, and 
compare it with the flux expected from a template star at the distance of the Orion Nebula. When an excess is measured, 
we attribute it to accretion. We then use our adaptation of the GuUbring et al. (1998) formula to derive the total accretion 
luminosity. The mass accretion rate can be estimated if both the radius and stellar mass are known. The stellar mass 
requires the assumption of an evolutionary model. In the following subsections we detail each individual step. 

4.1. Synthetic Photom,etry 

To calculate the magnitudes and colors of our template stars, we adopt the synthetic photometry package iraf/synphotJ 
In addition to the four HST filters, we consider Johnson V, Cousins I, and Bessell JHK filters. For most of the fol- 
lowing discussion we denote as UBVI filters those from ifST/WFPC2. We explicitely specify when referring to the 
Johnson- Cousins passbands. Synphot is also used to obtain zero magnitude fluxes from the spectrum of Vega, which 
is also part of the synphot library for spectral calibration. The zero-magnitude fluxes in our four HST UBVI fil- 
ters are, respectively: ZP(F336W)=3.32 x 10"^ ergs s"! cm'^ A^^; ZP(F439W)=6.65 x 10"^ ergs s"^ cm'^ A-^; 
ZP(F547M)=3.55 x lO'^ ergs s'^ cm'^ A"!; ZP(F791W)=1.18 x 10"^ ergs s^^ cm^^ 

The template spectra are taken from the Bruzual-Persson-Gunn-Stryker (BPGS) Spectral Atlas, which is an extension 
of the Gunn-Stryker (1983) optical atlas ranging from the UV to the near IR. The catalog contains several stars belonging 
to the young Praesepe and Hyades clusters that are representative of the stellar photosphere of our pre-main sequence 
objects. We select in particular the nine stars listed in Table 2. In order to have our templates with extinction Ay — 0, we 
derive their effective temperatures from the V — Ic colors using the relation provided by Hillenbrand (1997) in Appendix C. 
The spectral types, alos derived using the scale adopted by Hillenbrand (1997) in Table 1, closely match those found in 
the most recent literature, with the exception of BPGS59=BD-|-38 2457. This source has an only an historical spectral 
type K8, instead of the K4 value we found. In the absence of recent data, we attribute this discrepancy to an error in 
the original spectral classification. Since for M-type stars the color-temperature relations are quite uncertain, for stars 
later than MO we use the more recent scale of Reid and Hawley (2000). For the effective temperature vs. V — Ic color we 
derive the best fit relation 

logTofr = 816.31 x {V - ~ 8520.9 x {V - hf + 29629 x {V - h) - 34317. (1) 

^ WFPC2 calibration data available at www.noao.edu/staff/doIphin/wfpc2_calib/, updated on May 31, 2002 
^ Synphot manual available at www.stsci.edu/instruments/observatory/synphot.html 
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To locate the target stars in the HR diagram, we wih apply the bolometric correction to the V magnitude, BCy- In 
general we adopt the BCy vs. log Teg expression also provided in Appendix C of Hillenbrand (1997), except for M-type 
stars where we use the relation 

BCv = 0.27 - 0.604 x {V - 4) - 0.125 x {V - Icf (2) 

from Reid & Hawlcy (2000). 

To obtain the colors and bolometric corrections of intermediate spectral types wc interpolate between the template 
results using the effective temperature as a free parameter. Five target stars turn out to have effective temperatures lower 
than that of our coldest template, Gl 15B. We therefore extrapolate their colors using a second order fit to the colors of 
our four coldest templates. 

4.2. Red-leak in the F336W filter 

It is crucial to use a reliable model of the HST photometric system, especially for the F336W filter which is affected 
by red-leak (at 7500 A, it has 1% of the peak transmission at 3500 A). If the red-leak is ignored, it will lead to an 
overestimate of the observed U-band excess for stars with the reddest spectral energy distributions. This is shown in 
Figure 4, where we plot for the five reddest template stars the color index F336W-F439W in function of the extinction 
Av^ both for the Ry = 3.1 and Ry = 5.5 reddening laws (these two values are discussed in the next Section). While 
one would expect to have a redder color when the extinction increases, i.e. cmves rising to higher positive values, there 
is an anomalous "bluing" of the color index with the extinction (curves turn down). This behavior is well understood 
for WFPC2: according to the instrument handbook, the synphot predictions match the observations within 0.1% for 
all template stars. In general, however, the corrections to the photometry for the red leak produced with these models 
will depend on the accuracy of the temperature and extinction estimates. These uncertainties are dominated by stellar 
variability, and will be discussed in Section 5.1. Figure 4 indicates that the uncertainty on the reddening law can also 
play some significant role. 

4.3. Reddening 

The synphot task is also used to calculate the colors of the template stars in the presence of reddening. It is known that 
the reddening law in the direction of the Orion nebula is peculiar, with Ry = Av/E{B — V) ~ 5, instead of the usual 
Rv = 3.1 (Baade & Minkowski 1937, Johnson 1967). A higher ratio of total to selective extinction means that the dust 
extinction is relatively grayer, as expected for grains larger than those producing the standard interstellar reddening. It 
also means that a higher extinction is needed to produce the observed color excess. This effect increases the estimates of 
intrinsic stellar luminosity. 

To treat the two reddening laws in a homogeneous way, we use the Fitzpatrick (1999) prescription. In particular, for the 
standard interstellar reddening law we adopt the tabulated curve for Ry = 3.1. The reddening law for Orion is calculated 
using the IDL procedure FMRcurve.pro, also provided by Fitzpatrick (1999), with the parameters matching those in 
Table 5 of Fitzpartick & Massa (1990) for O^Ovi C (HD37022), i.e. Aq ^ = 4.635, 7 = 0.846, ci = 1.251, C2 = 0.033, and 
C4 = 0.186. The parameter R, which in the original Fitzpatrick's procedure is an input variable, is now calculated from 
the value of C2, providing R — 5.504. Once the reddening laws have been obtained, we attenuate the stellar templates by 
various amounts of extinction ranging between Ay = 0.1™ and Ay — 6™. To estimate the extinction properly, one must 
take into account also the possible contribution of the accretion spectrum, as discussed in the next subsection. 

4.4. Moss accretion and the Gullhring relation in the HST/WFPC2 F336W filter 
The relation of Gullbring et al. (1998) ("Gullbring relation" hereafter): 

log = imtlZ^ ■ log + 0.98tro?, (3) 

linking the Johnson U-band excess, Luj^cxccss to the total accretion luminosity, ^accretion, is based on observational data 
and is valid for the Johnson-U filter only. Our F336W filter is better suited to this type of study than the Johnson-U 
filter, as the F336W bandpass is entirely on the blue side of the Balmer jump, whereas the Johnson-U filter runs across 
it. This is an advantage, since, in general, the photosphcric emission of young stars drops or remains constant shortward 
of the Balmer discontinuity, while the ionized emission from the accretion material increases. In particular, whereas 
in the Johnson U-band the UV excess is higher for late spectral type stars than for early- type stars, due to partial 
cancelation of the emission and absorption from the ionized gas and hot photospheres, in our F336W filter the excess is 
expected to remain consistently high and increasing with the stellar temperature. The Gullbring relation was reproduced 
by Calvet & Gullbring (1998) on the basis of magnetospheric infall models (Uchida & Shibata 1984, Camenzind 1990, 
Konigl 1991). They assumed that the material, migrating from the disk to the stellar surface through accretion columns, 
produces shock emission that can be calculated with a 1-d model. In the wavelength range shortward of the Balmer 
discontinuity the spectrum is dominated by optically thin emission from the pre-shock gas and from the attenuated post- 
shock regions, whereas the Paschen and Bracket continua are mostly produced by the optically thick emission from the 
heated photosphere below the shock. The Gullbring relation predicted by this model has a slightly shallower slope than 
the observed one, with no evidence of dependence on the stellar spectral types of the underlying photosphere, at least in 
the range from K5 to M3. 
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Since our main goal is the calibration of the GuUbring relation to the F336W filter, we are mostly interested in the local 
behavior of the spectrum in the vicinity of the Balmer jump. Thus, we assume for simplicity that the excess is produced 
by a uniform and isothermal slab of ionized gas. Although a single slab cannot reproduce the entire accretion spectrum 
(Calvet & Gullbring 1998), it allows us to calculate a grid of models easily and to select among them those matching the 
original Gullbring relation within the errors. Our assumption is that the sub-set of models that reproduce the observed 
relation between L{7j,excess and -^accretion will also provide a reliable relation between LF336W,excess and ^accretion- In other 
words, we cross calibrate our slab models with the observational data in the region of the Balmer discontinuity. 

To generate the slab models, we use the photoionization code CLOUDY (Ferland 1996^) with a range of densities and 
temperatures matching those considered in previous studies (Batalha, Lopes, & Batalha 2001; Calvet & Gullbring 1998). 
In particular, the densities arc in the range logn = 14.6 — 15.2 and the temperatures in the range logT = 3.88 — 3.92. We 
consider the typical composition of HII regions as a baseline, but also explore values of 10 times higher or lower metallicity. 
The slab has a thickness of 10^'^ cm, and an excitation source of 10 Lq with a blackbody spectrum of Tgs = 10, 000 K, at 
a distance d = 10^^ cm from the slab. 

Each synthetic spectrum is convolved with both the Johnson-U filter and the F336W filter. The spectra matching the 
Gullbring relation within the errors are selected to generate the new relation between the accretion luminosity in the 
F336W filter and the total accretion luminosity. The revised Gullbring relation for the F336W filter is 

log (^h^^^'^ = (1.16 ± 0.04) • log j^ ^^33aM^,excess ^ ^ ^^ 24 ± 0.10). (4) 

The higher intercept of this relation (1.24 vs. 0.98) indicates that our F336W bandpass, narrower than the Johnson U 
bandpass, includes a smaller fraction of the accretion luminosity. On the other hand, the higher slope (1.16 vs. 1.09) 
indicates that the F336W bandpass is a more sensitive indicator of the accretion luminosity. 

The ionized gas responsible for the UV excess shortward of the Balmer discontinuity also affects the entire spectrum. 
Hartigan et al. (1991) have shown that the accretion continuum may contribute significantly to the flux in the V and 
R passbands, affecting the extinction estimates by ~ 50%. The emission in the Paschen continuum, in particular, may 
increase the measured flux in the I-band appreciably. To account for this contribution, we consider a representative slab 
model at T = 7600 K and repeat the procedure followed for the stellar templates, using synphot to predict the flux in 
all broad-band filters with differing amounts of reddening. Moving away from the Balmer discontinuity, our single slab 
models become less reliable, but still adequate to our purpose. 

4.5. Model fit and accretion luminosity 

To calculate the extinction, we use an iterative procedure. First, we determine Ay by comparing the observed vs. 
template colors. Then wc add our representative accretion spectrum, with the same amount of extinction just obtained 
for the pure photosphere, and match the observed U-band magnitude by adjusting a multiplicative constant that can 
be interpreted as a filling factor. A new extinction estimate is then obtained by comparing the observed data with the 
star -|- accretion spectrum, and the process repeated. With respect to the study of Hillenbrand (1997), which was based 
on the V-I color and assumed the standard Ry = 3.1 reddening law, we have the advantage of the extra B-band. A best 
fit to the BVI photometry provides a more reliable estimate of the extinction than the V-I color only, and possibly a hint 
on the best reddening law. On the other hand, the B-band data play in our case a special role. Since the B-band and 
the U-band images have been taken almost simultaneously, we expect the U-B color to be quite independent on source 
variability, being possibly affected only by flaring activity in the accretion process, but not by variability modulated by 
stellar rotation. In other words, to estimate the UV excess it is convenient to constrain the model to match exactly the 
observed B-band photometry. Our strategy is therefore to obtain the extinction from the best-fit to the observed BVI 
colors and the average stellar luminosity, and then normalize it to pass through the B-band point, providing in this way a 
most reliable estimate of the UV excess at the time of the observations. Technically, we use a least square estimator. To 
calculate the x^, we weight each deviation squared with the variances obtained by adding quadratically three contributions: 
the typical errors of the original measures in each band (reported in Table 1), the zero point errors quoted by Dolphin, 
and the standard deviations av{F547M) listed in the last column of Table 1, assumed to be representative of the stellar 
variability. In the four cases with only one available observation in the F547M filter, we assume that variability is present 
at the level of the mean standard deviation of the sample, ct = 0.211™. This procedure is carried out for both reddening 
laws. 

Once the extinction is known, the total stellar luminosity is obtained through the formula 

log (^^^ = O.4(M0,boi - Mtot.boi) = 0.4(4.75 -V- BCy + Ay + DM), (5) 

where we assume Mq ^hoi = 4.75'" for the solar absolute magnitude and a distance module to the Orion cluster, DM = 
8.36"*, corresponding to a distance of 470 pc. 
The accretion luminosity in the U-band, 

^U,excess = -^i/jobservation -^i/jexpected? 

refers to the luminosities corresponding to the dereddened magnitudes. The Gullbring relation for the HST U-band 
derived in the previous section (§4.4) provides the total accretion luminosity, ^accretion- We estimate the mass accretion 

® Updated on 2002 program and documentation for Beta 4 of Cloudy 96 available at http://nimbus.pa.uky.edu/cloudy/cloudy_96b4.html 
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rates by using the standard relation between M and iaccretiom 

TiV -^accretion-^^ /^\ 

^= 0.8GM. ' 

where the stellar radius, R^,, is obtained from the standard = AnRi^aT^ relation with cr — 5.67 x 10~^ ergcm"^ s^-'^ K^'*, 
and the stellar mass, is obtained from theoretical models of the pre-main-sequence evolution (§5.2). The factor of 0.8 
accounts for the fact that the material is assumed to be freely falling along the magnetosphere from a corotation radius 
Rmag ~ 5i?* (GuUbring et al. 1998). 



5. RESULTS 

5.1. Spectral Energy Distributions 

Figures 2 and 3 show the model fits to the observed magnitudes for the Ry = 3.1 and Ry = 5.5 reddening laws, 
respectively. For each star, we plot the measured magnitudes, the model spectrum (star+accretion) matching the BVI 
photometry (dotted line), the same model normalized to pass through the B-band photometry for the estimate of the 
U-band excess (solid line), and, when relevant, the accretion spectrum (dashed line). In Table 3 and Table 4 we list the 
corresponding results for the stellar luminosity (columns 3), the extinction (colmnn 4), the stellar radius (column 5), and 
the accretion luminosity (column 6), again for both reddening laws. The effective temperatures (column 2) have been 
taken from the Hillenbrand 1997 scale. 

On the highly compressed scale of Figures 2 and 3, the dotted lines nicely match the V- and I-band data. The 
displacement with respect to the solid lines, i.e. the deviation of the B-band with respect to the best fit, visualizes the 
range of variability between the average magnitude and the magnitude at epoch of our UV observations. In almost all 
cases the U-band flux appears strongly in excess with respect to the ideal reddened black-body curve, but this is largely 
due to the red-leak. The real accretion luminosity is represented by the distance between the U-band point and the solid 
line. When this is large, the standard accretion spectrum, normalized to account for the U-band flux, becomes evident at 
the bottom of the figure. The difference between the values of Ay derived from the two reddening laws is of the order of 
15%. The stellar luminosities are also affected in a similar way, due to Equation 5. 

To assess the reliability of our model fits we can refer to the parameter that is used to minimize the discrepancy 
between the model fit and the observed B, V, and I magnitudes, as described earlier. If one assumes that the deviations 
are approximately Gaussian, i.e. with zero mean and independent, then the distribution of the values should follow the 
theoretical reduced distribution, where v is the number of degrees of freedom. In our case, it is = 3 since for each star 
we typically have five bands (the WFPC2 BVI and the VI data from Hillenbrand (1997) minus two model parameters 
(reddening and normalization). Poor fits provide highly improbable x^ values, above some defined critical threshold. With 
40 sources, we can set the critical threshold at the 99% probability, X3(0.99) = 11.3. Above this value we should find 1% 
(i.e. none) of our stars, whereas we have approximately 50% of our sources, no matter which reddening law we adopt. 
According to the standard prescriptions, for these objects we should reject the hypothesis at the basis of the model fit, e.g. 
the quoted spectral types, or the adopted reddening law, or both. On the other hand, if systematic contributions to the 
noise are underestimated, a x^ behind the critical threshold may force wrong decisions. In our dataset, stellar variability 
is a prime suspect for loose fits, since our sparse sampling in time certainly underestimates the variability range. 

When the measured x^ > 11-3 and fT(F547M) < a = 0.211™, we make the assumption that variability has been 
underestimated. If we recalculate the x^ forcing in these cases (T(F547Ai") = ct, nearly half of the sources above the 
critical threshold can be recovered. The resulting x^ distributions, shown in Figure 5 together with the theorical xi 
distribution, show that the number of outliers above x^ — 6 (where the probability falls to approximately 10%) is still 
too large. Figure 5 clearly shows the limitations intrinsic to the dataset at our disposal. In particular, all sources with an 
anomalously high value of must be considered with great caution. For this reason, in the following we shall distinguish 
between a main sample, composed by sources with x^ < 11-3 that are expected to have a good fit, and a secondary 
sample with x^ > 11-3, having a less reliable fit. A Kolmogorov-Smirnov test on the two distributions of the observed x^ 
presented in Figure 5 indicates that the hypothesis that both are drawn from the same population must be accepted with 
a very high degree of confidence. In other words, the available data do not allow to select statistically one reddening law 
or the other. Since the reddening law affects the estimates of the mass accretion rate, we will continue to carry out our 
analysis for both laws. 

A last remark concerns the excess UV luminosity. Only in 7 cases {Ry = 3.1) or 9 cases {Ry = 5.5) there is no evidence 
of a UV excess. The lack of excess may be due to short term variability, related e.g. to flaring, which is not properly 
sampled in time-scale of our UB observations. In Figure 6 we plot the histogram of the U-band excess luminosity for 
the two different reddening laws. The strong asymmetry to the positive side, even at the lowest signal levels, indicates 
that the F336W excess is real. Figure 2 indicates that only a systematic underestimate of the extinction could artificially 
produce an excess emission in the UV, especially for the most reddened stars. Given that the main uncertainty in our 
best fit is caused by random stellar variability, we tend to rule out major systematic errors. Since the ratio between the 
accretion luminosity and the total stellar luminosity, corresponding to the veiling factor of the stellar absorption lines, is 
generally small, the positions of the stars in the HR diagram remain almost unaffected by the removal of the accretion 
contribution. 



HST/WFPC2 Study of the Trapezium Cluster 



7 



5.2. The HR diagram 

The location of our 40 stars in the HR diagram is shown in Figure 7. The four plots have been obtained assuming the 
Rv = 3.1 (loft column) and Ry = 5.5 (right column) reddening laws, and two different sets of theoretical models. On the 
top row we use the latest version of the isochrones and evolutionary tracks of D'Antona & Mazzitelli.^ (1998, hereafter 
DM98), whereas on the bottom row we use the results of Palla & Stahler (1999, hereafter PSF99) The different symbols 
identify the sources according to their reliability class: filled dots are used for the more reliable sources of the main sample, 
small dots with circles are used for the remaining more uncertain sources. On each figure, the two diagrams located at 
the bottom at log(Toff ) = 3.694 (K2 spectral type) and at the top at log(Teff) = 3.502 (M3.5 spectral type) represent how 
different types of errors change the position of a stars in the diagram. Namely, 1) the vertical bar represents the ±lcr 
error associated to the average photometric uncertainty AV = 0.211'", that we attributed to stellar variability; 2) the 
horizontal bar represents the error on the effective temperature, reported by Hillenbrand (1997) to be equal to ± 1 subclass 
(corresponding to A log Teg « ±0.011) for stars K7 (i.e. logTeff = 3.602)' and later, and ±1/2 class (A log Teg w ±0.05) 
for earlier stars; 3) the diagonal bar represents the error in the bolometric correction that is introduced when a inaccurate 
spectral type is assumed. Since for our late spectral types the bolometric correction increases (becomes more negative) 
with the spectral type, a star of a given magnitude will be more luminous if a lower temperature is attributed to it. A 
lower temperature will also make the spectral type redder, thus reducing the estimated extinction; this compensates for 
the error on the bolometric correction by an amount that depends on the extinction. In the plot we refer to the case 
without extinction. 

Using theoretical evolutionary models, one can interpolate between evolutionary tracks and isochrones to derive the 
mass and the age of each individual star. In our case, we are mostly interested in the mass value as it enters in Equation 6 
for the estimate of the mass accretion rate."'^'' The results, also listed in Table 3 and 4, are strongly model dependent 
for a variety of reasons as discussed, e.g., by Baraffe et al. (2002). The two particular models we have considered are 
characterized by different choices in the initial conditions (stellar radius), but in addition other physical processes, like 
chemical composition, convection efficiency, surface gravity, and even the accretion rates that we are trying to measure, 
contribute to provide discrepant results. According to Baraffe et al. (2002), at present there is no reliable theoretical 
model for ages < 10^ yr. Even if absolute ages and masses cannot be unambiguously determined, the comparison of the 
results obtained using different models may provide useful insights. 

Using DM98 models, masses range from to 4.2 M© to less than 0.08 Mq, the H-burning limit. The values, however, 
may vary by a factor of ~ 2.0 if one uses different PMS models, especially in the range ~ 0.6 — 1 Mq. For example, 
source 9209 has logT = 3.58 and log{L/ Lq) — 3.078 in the reference b) case. The corresponding mass is M = 0.26 Mq 
using D'Antona and Mazzitelli tracks, or M = 0.57 Mq using PS99 tracks. The stellar masses of the sources above the 
0.1 Myr isochrone (for DM98) or above the birthline (for PS99) have been estimated using the mass of a source lying on 
the 0.1 Myr isochrone or birthline at the same effective temperature. On the other hand, there are five sources that lie 
below the 0.1 Mq evolutionary track of PS99. These most interesting brown dwarf candidates fall, in the PS99 case, in 
our secondary sample. Given the uncertainties, we did not attempt to extrapolate our fitting process behind the 0.1 Mq 
track. With these provisions, the average mass of our sample turns out to be 

1. D'Antona & MazziteUi: < log M(M0) >= -0.52 ±0.35 

2. Palla & Stahler: < log M(M0) >= -0.371 ± 0.37 

essentially regardless of reddening law. 

Concerning ages, the evolutionary tracks and isochrones of D'Antona & Mazzitelli have been used by Hillenbrand (1997) 
to derive a cluster age less than 1 x 10^ yr for masses between 0.01 and 2.5 M0. Using near-IR data and updated D'Antona 
& Mazzitelli (1998) models, Hillenbrand & Carpenter (2000) found that a star formation rate constant between 3 x 10'' 
and 3 x 10^ years in logarithm of the age bins closely matches the distribution of stellar ages. Luhman et al. (2000) 
compared their HST/NICMOS near infrared observations and ground-based K-band spectra with D'Antona & Mazzitelli 
(1998) models and foimd a median age of 4 x 10"'' years. PS99, on the other hand, used their models on the same dataset of 
Hillenbrand (1997), deriving a best fitting single age of 2 x 10^ yr, with star formation starting at low level some 10^ years 
ago and accelerating to the present epoch (but see Hartmann 2001 for a discussion of the systematic errors). PS99 models 
use deuterium burning to control the protostellar radius during the initial accretion phase. This sets a well defined locus 
(birthline) in the HR diagram. The birthline shown in Figure 7 has been calculated for an assumed mass accretion rate 
of M = 1 X lO-'"^ Mgyr-i. 

Considering the logarithm of the ages, as recommended by Hartmann (2001), we obtain ages compatible with previous 
estimates based on these families of models: 

1. DM98+iiy=3.1: logi(yr) = 5.76 ±0.90 (average); logt(yr) = 6.02 (median) 

2. DM98+i?y=5.5: logi(yr) = 5.51 ± 1.07 (average); logi(yr) = 5.79 (median) 

3. PS98+i?y=3.1: logt(yr) = 5.99 ± 0.55 (average); logt(yr) = 6.07 (median) 

^ available at http://www.mporzio.astro.it/ dantona/ 

In Section 6 we will compare the mass accretion rates in Orion with those in Taurus, estimated by Hartmann et al. (1998) using the DM98 
family of models. 
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4. PS98+i?v=5.5: logt(yr) = 5.80 ± 0.58 (average); logt(yr) = 5.94 (median) 

If the best assumptions are those that minimize the scatter with respect to an average cluster's isochrone (Hartmann 
2001), then PS99 models seem preferable. PS99 models also predict isochrones that apparently show a better match to 
the overall distribution of stars in the HR diagram. In particular, in the Ry = 5.5 case the maximum age is consistently 
close to the 3 Myr isochrone, whereas the D'Antona & Mazzitelli isochrones show a correlation between ages and mass, 
with more massive stars appearing younger. This trend was already noticed by Luhman et al. (2000), Hillenbrand & 
Carpenter (2000) and Muench et al. (2002). The age distribution in this case (Figure 8) shows some evidence of the 
elusive "pileup" of stars at the maximum age expected when star formation is concentrated in a relatively short episode 
(Hartmann 2001). The increase of the number density of stars is roughly linear with the logarithm of the age, and can 
be approximated by the relation logdA^ = 3(logi — 4). The observed logarithmic decline actually represents an increase 
of the star formation rate with time, similar to the "acceleration" suggested by PS99. In this scenario, the Orion cluster 
can be considered as an active site of star formation. 

5.3. Mass Accretion Rates 

Following the method described in Section 4.5, the mass accretion rate, M, can be estimated for approximately 3/4 
of the stars. The values, listed in Table 3 and 4 for both evolutionary models and reddening laws, range between 
~ 2 X 10^^ — 1O~'^^M0 yr~^. For the remaining 1/4 of stars, Al cannot be obtained either because there is no mass 
estimate (the star lies in the HR diagram out of the model range) or because there is a deficit of UV flux with respect to 
the stellar photosphere. Stars of the main sample have typically mass accretion rates lower than those in the secondary 
sample. Moreover, stars with UV deficit are typically part of the main sample. In conclusion, when a star is well measured, 
the mass accretion rate is always low, or absent. The average and median values for the main sample are: 

1. DM98 + Rv ^ 0.3: avcrage=3.59 x 10"^ Mq yr'^; median=7.4 x 10"!° Mq yr"! 

2. DM98 + Rv ^ 5.5: average=3.12 x 10"^ Mq yr'^; median=1.0 x 10"" Mq yr"! 

3. PS99 + Rv ^ 3.1: average=3.90 x 10"" Mq yr'^; median=6.0 x 10"^° Mq yr'^ 

4. PS99 + Rv = 5.5: average=2.39 x 10"" Mq yr"i; median=8.6 x 10"^° Mq yr"! 

In Figure 9 we show the same four HR diagrams of Figure 7, using for each star a circle with diameter proportional to 
the logarithm of M, if measured. Thick and thin circles represent main and secondary sources, respectively. Figure 9 
shows no obvious correlation between mass accretion rates and the position of the star in the HR diagram. A direct 
comparison of the four plots reveals how different reddening laws and PMS evolutionary models affect the estimates of M. 
The uncertainty on the reddening law affects the measure of M through the U-band extinction, the red-leak correction 
and the stellar radius, derived from the stellar luminosity. Evolutionary models provide different values of the stellar 
mass, that directly enters in the estimate of M through Equation 6. Table 3 and 4 show that the average scatter between 
different estimates of accretion rates is of the order of 30%. 

Besides the reddening laws and theoretical models, there are other sources of uncertainty. Since the Orion cluster is 
relatively far and compact, the relative distance of cluster members has a negligible effect on the scatter of the luminosities 
and mass accretion rates, assuming the cluster depth of the order of the projected cluster size, ~ 1 pc. A major 
systematic uncertainty is the average distance of the Cluster, which is known to within a ^ 20% error. It ranges from 

~ 400 - 480 (Warren & Hcsser 1977) to 470 ± 80 (Genzel et al. 1981). Changing the distance from 470 pc to 430 pc 
decreases the stellar luminosities by a factor « 20% and the mass accretion rates by a factor « 28%, on average. Another 
systematic error is introduced by the assumed value for the radius at which magnctospheric infall begins, i?mag- The value 
we assumed, i?,nag = 5, is equal to that assumed by Hartmann et al. (1998) for their study of the Taurus and Ophiucus 
region, but other choices are possible. MuzeroUe et al. (2003) have recently assumed for instance i?niag = 3, on the basis of 
their Ha emission line models. This choice would lead to mass accretion rates lower by ^ 20%. Concerning photometric 
errors, they also contribute to the uncertainties on the accretion rates. An errors of ~ 0.05"* in the U band magnitude 
adds ~ 5% of uncertainly. In general, measurement errors provide a negligible uncertainty with respect to systematic 
effects we have discussed. On average, our results should be correct to within a factor of 2. 

6. DISCUSSION 

6.1. UV excess as a tracer of mass accretion 

We have selected a sample of young {t < 3 Myr), low mass {M ~ 0.1 — 1 Mq) stars with evidence of mass accretion 
(M ~ 10~^ — 10~^^ M0yr~^) in ~ 75% of the cases. The interpretation of UV excess as entirely due to mass accretion rates 
can be questioned. Disk photoevaporation creates envelopes surrounded by a photoionization front facing the ionizing 
stars (see later). The recombination flux emitted by these low density ionization fronts should be quite different from that 
produced in the shock excited high-density accretion column, but our broad-band photometry does not allow to distinguish 
between these two processes. To mitigate possible contamination, we have rejected from our sample all sources appearing 
clearly extended on a visual inspection. A second check is offered by Figure 10, where the location of our sources in the 
Trapezium core is shown. Each star is represented by a circle with diameter proportional to the excess UV luminosity. If 
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radiation by Trapezium stars was to be the main source of stcUar envelope ionization, one would expect to find a strong 
gradient of the UV flux with the distance from the core of the Nebula, which is not observed. Within the area covered 
by our study there is no clear trend between the excess UV luminosity and the projected distance from the center of the 
nebula. In any case, any contamination from the envelope would increase the UV excess, forcing us to consider the values 
wc have found as upper limits. Since the accretion rates appear already anomalously low, this effect would just strengthen 
our conclusions. 

The interpretation of the UV excess as actually due to the accretion process is strengthened by the strong correlation 
with the near IR excess, represented in Figure 11 through the K-band excess. This relation has been quite elusive in 
the past, but in our case turns out to be evident. This because we can rely on accurate photospheric subtraction, rather 
than on color indexes. There scatter is still large, and category 1 stars with reliable spectral energy distributions show an 
excess stronger in the near-IR than in the UV. One must take into account that the UV and near IR emission come from 
different regions, and at small distances from the star the disk orientation plays a dominant role shadowing the inner disk 
region. A more detailed study of the IR excess emission may lead to constrain the inner holes and reduce the uncertainties 
on the assumed i?mag- In any case, the presence of siibstantial IR excess indicates that possible inner holes have small 
size, coherently with a scenario where the accretion process proceeds across the disk down to the magnetospheric radius. 
It is interesting to notice that if the growth of protoplanets is traced by the development of large disk gaps, our disks 
appear to be in an earlier evolutionary phase. 

6.2. Low mass accretion rates 

In Figure 12 we plot the mass accretion rates derived assuming the DM98 and PS99 models and Ry = 3.1 and Ry = 5.5 

reddening laws, versus the stellar ages. Black filled and open circles are \iscd for main and secondary sources, respectively, 
and their size is proportional to the stellar mass. To reduce the uncertainties related to the definition of the zero age, we 
have considered only stars with ages older than 0.1 Myr. 

The red circles refer to the Taurus sources studied by Hartmann et al. (1998) assuming Ry = 3.1 and the DM98 tracks, 
and using, for most sources, the same UV excess technique adopted in our study. For their sample Hartmann et al. (1998) 
found a decrease of the mass accretion rate with stellar age. In our case, the possibility of a correlation depends on the 
assumptions. A Spearman rank-order test, which being non parametric is independent on the actual values of age and 
accretion rate, shows that chance of a random correlation varies in the four cases according to: 

1. DM98 + Rv = 0.3: P=3.7% (main sample), P=3.6% (all stars) 

2. DM98 + Rv = 5.5: P=0.9% (main sample), P=60% (all stars) 

3. PS99 + Rv =3.1: P=0.02% (main sample), P=0.09% (all stars) 

4. PS99 + Rv = 5.5: P=13% (main sample), P=18% (all stars) 

If we restrict to the more reliable sources of the main sample, evidence of correlation seems especially strong in the case 
of PS99 + Rv = 3.1. In this case, a power law fit to the entire sample shows that M decreases as M(t) = Kt^^, with 
rj = 1. The scatter, however, is large, with mass accretion rates that differ by about two order of magnitudes at the same 
stellar age. For comparison, Hartmann et al. (1998) find 1.5 < ry < 2.8 for Taurus and Cha I associations, and comparably 
large scatter. A more quantitative analysis of the correlation is complicated by the fact that the age and accretion rate 
are intrinsically correlated, being both derived from the stellar luminosity. 

The steeper decline with time (higher 77) for the Taurus sample is not the only difference with our Orion results. 
Figure 12 clearly shows that mass accretion rates in Orion are systematically lower: in the previous section we have 
quoted average values ~ 10~^ Mq yr~^, whereas Hartmann et al. (1998) find for Taurus an average of 10~^ Mq yr~^. 
Since Hartmann et al. (1998) sample contains a similar number of stars with ages comparable to our sources, we tend 
to reject the hypothesis that all of our stars have been observed by chance in a quiescent phase. The same mean value 
of Hartmann et al. (1998) was also found by RebuU et al. (2000) in their study of the mass accretion rates in the Orion 
Nebula flanking flelds, four 45' x 45' flelds centered ~ 0?5 east, west, north, and south of the core area covered by our 
study. Even if we consider the stars younger than 0.1 Myr, not plotted in Figure 12, the mass accretion rates remain of the 
order of 10~^ Mq yr"-"^. The flatter decay curve we find in Orion therefore appears related to a recent systematic leveling 
off of the mass accretion rates, while one would have expected finding some significant increase at the young ages probed 
by our study. PS99 have shown that the locus of the Trapezimn stars in the HR diagram is nicely enveloped, at the high 
luminosity side, by the birthline calculated for a protostar with accretion rate of 1O~^M0 yr~^ and deuterium-to-hydrogen 
ratio D/H=2.5 x 10~^. The locus of the birthline in the HR diagram does not depend upon the initial mass and radius 
of the collapsing protostellar core, but depends on the mass accretion rate: given two cores with the same mass, the one 
having lower mass accretion rate will be more compact. Due to the higher density, it will ignite and deplete its deuterium 
content earlier and be less luminous (Stabler 1988). The M ~ 1O~^M0 yr~^ value of PS99 is compatible with stellar 
collapse studies (Masunaga, Miyama, & Inutsuka 1998) and provides a good agreement with the observations, not only 
in what concerns the birthline, but also the zero-age main-sequence locus and the evolutionary tracks (Palla & Stabler 
2001). 

Stars on the 10~^ Mq yr~^ birthline may show lower mass accretion rates. Hartmann, Cassen, & Kenyon (1997) have 
found that if the accretion occurs mainly from an equatorial disk (rather than from a spherical envelope as in PS99), 
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a sudden drop in the mass accretion rate has no immediate influence on the photospheric effective temperature and 
luminosity. The residual time spent on the birthline, or in its vicinity, can be easily estimated, since the energy available 
from Deuterium fusion in a core of mass, Mg, is Lt = 1 Lq{1.5 x lO^yr) (Mc/Mq). A completely convective, non accreting, 
Mc = O.2M0 PMS star on the 10~^ Mq yr~^ birthline will need w 10® yr to resume contraction. This time drops by 
about an order of magnitude at Mc ~ 0.5 Mq, becoming negligible at higher masses. Therefore, in order to maintain in 
the vicinity of the birthline a group of stars with mass spread larger than an order of magnitude and very small mass 
accretion rates, one must assume that the main accretion phase has been recently terminated for all sources, independently 
on their ages. In principle, this can be simply due to the spontaneous exhaustion of disk material available for accretion. 
However, it seems rather artificial to have all disks depleted almost at the same time across the cluster. A more attractive 
possibility seems to attribute the drop of mass accretion to a common trigger event caused by an external agcnit. 

It is known that disks in Orion are photo-ablated by the UV radiation of OB stars, the "proplyd" phenomenon first 
found by Churchwell et al. (1987) and later confirmed by spectacular HST images (O'Dell, Wen, & Hu 1993, O'Dell 
& Wen 1994). Churchwell et al. (1987) estimated mass-loss rates of about 10~^ Mq yr~^, whereas Henney & O'Dell 
(1999) more recently have measured mass-loss rates as high as 10~^Mq yr~^. Thus, when a protostar is exposed to the 
ionizing flux of a new-born OB star, the disk mass decreases rapidly with time. The similarity solution of Hartmann et al. 
(1998) provides a theoretical link between the disk mass and the mass accretion rate. In Figure 12 we show the M-age 
relations expected for viscous disk accretion. The solid line corresponds to the fiducial model with r] = 1.5, stellar mass 
Mi, = 0.5 M0, radial disk scale factor i?i = 10 AU, outer disk temperature Td2 = 10 K, viscosity parameter a = 0.01, 
and initial disk mass Md,{0) = 0.1 Mq. The dashed line refers has the same parameters except for = 0.1 M©, whereas 
the dotted line has M^, = 0.1 Mq and Md(0) = 0.01 Mq. In general, the mass accretion rate is proportional to the disk 
mass and decreases with time. It is important to note that this model refers to the disk mass accretion, Md, whereas 
our observations refer to the stellar mass accretion, M^,. Mass conservation implies that part of the accreted disk mass 
is lost through the outflow activity at a rate My,. The standard model for magneto centrifugally driven winds predicts 
Md/Mu, « 3.5 (Shu et al. 1994), thus Md = 1.4M*, i.e. disk and stellar mass accretion rates are of the same order. In 
Orion, mass loss is typically traced by "microjets" with modest mass-loss rates, around lO^^M© yr~^ (Bally, O'Dell, & 
McCaughrean 2000). 

Since the mass loss from the disk evaporation dominates all other mass transfer mechanisms, it is tempting to speculate 
that the reduction of disk mass resulting from the exposure of a circumstellar disk to the UV radiation of newly born 

OB stars regulates the mass accretion rate through the disk, and therefore to the star. The similarity solution, based on 
assumptions like thin disk geometry and negligible external heating, suggests a possible link: the mass accretion rate is 
proportional to the disk mass. The similarity solution, however, can hardly be directly applicable to real photoevaporated 
disks. In a previous paper we have shown that when the circumstellar disk becomes exposed to the UV radiation, the disk 
flaring angle (i.e., thickness) and temperature show a dramatic increase, depending mostly on the distance and orientation 
of the disk with respect to the main heating source (Robberto, Beckwith & Panagia 2003). An increase of the outer disk 
temperature raises the efliciency of internal energy dissipation, supporting higher mass accretion rates. Further theoretical 
investigations are required to understand the interplay between disk photoevaporation and mass accretion. 

6.3. Implications for the star formation history and IMF 

If the mass accretion rates are affected by the disk photoevaporation phenomenon, one may expect to find systematic 
differences between the stellar population in the core of the Orion Nebula and other regions, like Taurus or the outskirts 
of the Orion Nebula itself, where the ionizing flux is negligible. For PMS sources still in the disk accretion phase, the 
sudden decrease of the mass accretion rate causes an premature end of the stellar mass build-up. Whereas the final stellar 
mass is fixed in the early stages of protostar formation, through competitive collapse/fragmentation phenomena regulated 
by magnetic fields or supersonic turbulence, most of the stellar material is accreted from the circumstellar disks at a 
rate which decreases with time. There is growing evidence that low mass stars and sub-stellar mass objects share this 
type of evolution (Muzerolle et al. 2003), possibly on longer time scales. These may especially be affected by the sudden 
disk evaporation, growing to a final mass lower than that they would have attained if the star formation process would 
have proceeded in a more quiet environment. Low mass objects, therefore, may remain "dwarfed" by the sudden disk 
dissipation, resulting in a relative overabundance of low mass stars and brown dwarfs. Various authors, and in particular 
Luhmann et al. (2000), have recently reported an overabundance by a factor of 2 in brown dwarfs in Orion relative to 
Taurus. At the same time, the overabundance of low mass stars and brown dwarfs should be compensated by a depletion 
of intermediate mass stars. In this mass range, comparison are affected by small number statistics. Still, it is intriguing 
to notice that Hillenbrand (1997) foimd a flattening of the IMF in the core of the Trapezium cluster at masses lower 
the 0.6 Mq, whereas the overall Orion stellar population follows the Salpeter law down to the completeness limit of her 
survey, 0.1 Mq. This flattening could be the final outcome of the disk depletion, causing all stars cascading in lower mass 
bins. The ultimate consequence of this scenario is that the low-mass end of the initial mass function could be modulated 
by the recent star formation history, namely by the formation of the OB stars in the cluster. 

7. CONCLUSION 

We have performed the first quantitative analysis of the mass accretion rates in the core of the Trapezium cluster, 
focusing our attention on 40 unresolved stellar sources with known spectral types and accurate WFPC2 UBVI-band 
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photometry. We have estimated the ultraviolet excess considering two different reddening laws, and derived the stellar 
parameters and mass accretion rates using the D'Antona & Mazzitelli (1998) and Palla & Stahler 1999) computations. 
Approximately 75% of the sources show appreciable excess luminosity in the U-band, that we attribute to accretion. 
We used the known correlation between the U-band excess and the total accretion luminosity (GuUbring et al. 1998), 
recalibrated to our photometric system, to estimate the mass accretion rates, all found to be in the range 10~* — 
1O~^^M0 yr~^. We find some evidence for a decrease of M with the stellar age for sources older than 0.1 Myr, but 
in general the mass accretion rates appear to be lower than those measured in Taurus or in the flanking fields of the 
Orion Nebula. We suggest that the photoionization generated by the Trapezium OB stars causes a drop of the disk 
mass accretion rate. Low mass stars therefore conclude their pre-main sequence evolution with lower masses, and the 
Initial Mass Function turns out to be affected by the rapid evolution of the most massive cluster's stars, with a surplus 
of "accretion aborted" stars or brown dwarfs, and a deficit of intermediate mass stars. This trend is in agreement with 
recent observations of the IMF in the Trapezium cluster. 

The authors are indebted to Ed Fitzpatrick for discussions on the reddening law, to Francesco Palla and Keivan Stassun 
for their comments on an earlier version of the manuscript, to Charlie Lada and Michael Meyer for discussions, and to 
B. Hilbert and M. Richardson for their collaboration. J. S. and G. M. C. have been supported by the Summer Student 
Program of the Space Telescope Science Institute. Support for J. S. has been also provided by HST/DDRF grant 82316. 
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Fig. 1. — Scatter of the measured F547M (V-band) magnitudes vs. the mean values. Filled dots refer to WFPC2 data, crosses refer to the 
V-band data of Hillenbrand (1997). 
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2 r 




Fig. 4. — Effect of the red-leak in the F336W filter, shown has a variation of the F336W-F439W color vs. extinction Ay for our template 
stars. Color is bluer to the bottom. Solid lines are used for Ry = 5.5, dotted lines for R = 3.1. Labels at the right edge refer to the logarithm 
of the stellcir effective temperature. 




Fig. 5. — Histogram of the minimumx^ distribution, for both reddening laws. The solid lines refer to the theorical X2 reduced to 3 degrees 
of freedom. 
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Fig. 6. — Histogram of the F336W excess, after subtraction of the photospheric continuum. Bluer colors are to the right. Open areas are 
used for the entire sample of 40 stars, hatched areas for sources of the main sample. 



R,= 3.1 R,=5.5 




3.80 3.70 3.60 3.50 3.40 3.70 3.60 3.50 3.40 

LogT.,, LogT.,, 

Fig. 7. — HR diagrams obtained for different reddening laws (Rv = 3.1, left colum; Ry = 5.5, right column) and evolutionary models 
(DM98, top row; PS99, bottom row). Filled circles are used for category 1 (most reliable fit) sources, open circles with inner dot arc used for 
category 2 sources (less reliable fit). Ages and masses arc indicated with conventional notation, with "BL" in the PS99 diagrams indicating 
the 10~^Mq yr~^ birthline. The diagrams at logToff = 3.70, log(L/LQ) = —2.8 and logT^ff = 3.50, log(L/L0) = 1.5 represent the typical 
errors resulting from spectral classification, stellar variability, and bolometric correction (see text). 
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Fig. 8. — Histogram of the age distribution. Open areas axe used for the entire sample of 40 stars, hatched areas for sources of the main 
sample. Top row: DM98 model; bottom row: PS99 model; left column: Ry = 3.1 reddening law; right column: Ry = 5.5 reddening law. 
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Fig. 9. — HR diagrams with mass accretion rates. This figure, similar to Figure 7, represents each star with an open circles with diameter 
proportional to the logarithm of the mass accretion rates. Asterisks are used for sources with UV deficit. 
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Fig. 10. — Location of the accreting protostaxs in the Trapezium. Open circles, used for category 1 sources, are proportional to the logarithm 
of the F336W excess luminosity, estimated assuming Ry =3.1 reddening law; plus (+) signs are used for category 2 sources; times (x) signs 
are used for sources with deficit of F336W emission. Filled dots represent cluster's OB stars, with the five Trapezium stars clearly visible at 
the center of the figure. 
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Fig. 11. — Infrared (K-band) vs. UV excess, calculated assuming the Ry = 3.1 reddening law. Filled circles indicate category 1 sources, 
times (x) signs indicate category 2 sources. The error bars on the UV excess represent the uncertainty associated to the U-band photometric 
errors. 
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Fig. 12. — Mass accretion rate vs. age for the main sample. The four diagrams represent different choices for the reddening law {Rv = 3.1, 
left column; Ry = 5.5, right column) and for the evolutionary models (DM98, top row; PS99, bottom row). Dark Glled circles: main sources; 
open circles: secondary sources; red filled circles mass accretion rates estimated by Hartmann et al. (1998) for the Taurus cluster. Solid 
line: similarity solution from the fiducial model of Hartmann et al. (1998); Dotted line: similarity solution, with = 0.1 Mq; Dashed line: 
similarity solution, with M* = 0.1 M0 and Ma{0) = 0.01 Mq. 
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Table 2 

Template stars from the Bruzual-Persson-Gunn-Stryker Spectral Atlas 
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Name 
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Sp.T.'' 
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'^spectral types and effective temperatures derived 
from the V — Ic color according to Hillenbrand (1997) 
(see text). 



HST/WFPC2 Study of the Trapezium Cluster 



23 



Table 3 

Stellar and Mass Accretion Parameters for Rv = 3.1 



Palla fc Stabler (1999) D'Antona fc Mazzitelli (1998) 
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Table 4 

Stellar and Mass Accretion Parameters for Rv = 5.5 



Palla fc Stabler (1999) D'Antona fc Mazzitelli (1998) 
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